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Dynamical quantum phase transitions (DQPTs) at critical times appear as non-analyticities dur¬ 
ing nonequilibrium quantum real-time evolution. Although there is evidence for a close relationship 
between DQPTs and equilibrium phase transitions, a major challenge is still to connect to fundamen¬ 
tal concepts such as scaling and universality. In this work, renormalization group transformations 
in complex parameter space are formulated for quantum quenches in Ising models showing that 
the DQPTs are critical points associated with unstable fixed points of equilibrium Ising models. 

Therefore, these DQPTs obey scaling and universality. On the basis of numerical simulations, sig¬ 
natures of these DQPTs in the dynamical buildup of spin correlations are found with an associated 
power-law scaling determined solely by the fixed point’s universality class. An outlook is given on 
how to explore this dynamical scaling experimentally in systems of trapped ions. 

PACS numbers: 05.70.Ln,64.60.Ht,73.22.Gk 


Introduction In equilibrium phase transitions are of 
fundamental importance both for the theoretical under¬ 
standing of physical systems as well as for applications. 
In this context, continuous phase transitions are of par¬ 
ticular interest because they exhibit scaling and univer¬ 
sality |Ti. These fundamental concepts are intimately 
connected to renormalization group (RG) theory and 
the associated fixed points. Dynamical quantum phase 
transitions (DQPTs) during quantum real-time evolution 
have emerged as a nonequilibrium analogue to equilib¬ 
rium phase transitions where Loschmidt amplitudes 

g(t) = (i/jo\e~ iHt \ii , o) ( 1 ) 

become nonanalytic as a function of time m- Here, \ip 0 ) 
is an initial pure state and H the Hamiltonian driving 
the coherent time evolution. By now, DQPTs have been 
discovered in a variety of different contexts EKES and 
indications for a close relationship between DQPTs and 
equilibrium phase transitions have been found EUMD- 
0®. But still, a major challenge is to connect to funda¬ 
mental concepts such as scaling and universality. 

In this work it is shown for the first time that 
DQPTs obey scaling and universality. For that purpose, 
Loschmidt amplitudes for quantum quenches in one- and 
two-dimensional Ising models are mapped exactly onto 
equilibrium partition functions at complex couplings for 
which RG transformations in complex parameter space 
are formulated. As a main result, DQPTs are criti¬ 
cal points on the attracting manifold of the unstable 
fixed points of this RG. Therefore, Loschmidt amplitudes 
satisfy a scaling form with exponents determined solely 
by the underlying universality class. Moreover, numeri¬ 
cal evidence is provided that the critical phenomena in 
Loschmidt amplitudes are related to dynamical power- 
law scaling in spin-spin correlations. An outlook is given 
on how to verify this scaling experimentally in systems 



t/t c 

Figure 1: (color online) Dynamical quantum phase transitions 
in the Loschmidt echo rate function A(t) = — N^ 1 log[|C?(t)| 2 ] 
after quenches in the ID Ising chain. The nonanalytic kink 
structure of A (t) is a direct consequence of the universal ex¬ 
ponents of the underlying fixed point, see Eq. §. 

of trapped ions within current technology. 

Universality and scaling of DQPTs will be studied for 
quantum quenches in transverse-field Ising models: 

L 

H(h) = - J2 -h^f- (2) 

(Im) l—l 

with er“, a = x,z, Pauli matrices on lattice site l = 
1, ...,L and L the total number of spins. While in 
one dimension (ID) the nearest-neighbor (NN) coupling 
Jim > 0 is taken as uniform Ji m = J, in two dimensions 
(2D) an anisotropic square lattice is considered with cou¬ 
plings J within the rows and J± along the columns. 

The Ising model supports DQPTs both in ID BUTO t m ] 
and also in 2D [T2|. Fig. [l] shows DQPTs for quantum 
quenches in the ID case. Throughout this letter, plots 
of the DQPTs will be given in terms of the Loschmidt 
echo C(t) = \Q(t)\ 2 quantifying the magnitude of Git). 
Specifically, due to the large-deviation scaling of C(t) [2] 
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mm it is suitable to introduce its rate function 

= (3) 

which, in contrast to C(t ), is intensive in the thermody¬ 
namic limit. As in Fig. [lj quantum quenches from fully 
polarized initial states 

l^) = 0H)H^ = ^[lt)i + U)j] (4) 


will be considered, i.e., ground states of initial Hamilto¬ 
nians Hq = H(h oo). 

In the following, it is the aim of the present work to 
relate the DQPTs and therefore the nonanalytic struc¬ 
ture of A (t) to scaling in the vicinity of unstable fixed 
points. For that purpose, it will be shown that Losclimidt 
amplitudes in Eq. 0 can be mapped onto equilibrium 
partition functions of classical Ising models at complex 
couplings for the considered parameter regime. In or¬ 
der to address scaling and universality, real-space dec¬ 
imation RGs, exact in ID and approximate in 2D, are 
formulated. The most important result of the analysis of 
the RG equations is that DQPTs are critical points flow¬ 
ing to unstable fixed points of equilibrium Ising models 
implying universality and scaling. Notice that singular 
behavior can also occur in the Fourier transform of the 
Loschmidt amplitude [15HT8] which, however, is of differ¬ 
ent nature. 

Equilibrium partition functions.- Let us first con¬ 
sider vanishing final transverse fields with a Hamiltonian 
H(h = 0) = — a i a m • Then, as will be outlined 

below, the Loschmidt amplitude Q{t) becomes: 


m = ^ 


itT, 


(i m ) 


JlmCTl 


(5) 


with Tr denoting the trace over Hilbert space. Remark¬ 
ably, this is nothing but the equilibrium partition func¬ 
tion of a classical Ising model at complex coupling iJ , 
inverse temperature t, and g(t) = — iV _1 log[C/(i)] the as¬ 
sociated free energy density (apart from an overall tem¬ 
perature normalization). The above identification follows 
by using two properties. First, \ifo) in Eq. Q can be 
written as an equally weighted superposition of all spin 
configurations in the erf basis. Secondly, H{h = 0) does 
not induce spin flips and therefore only the diagonal ma¬ 
trix elements contribute which is nothing but the trace. 
In analogy to the equilibrium case it will be suitable to 
introduce dimensionless couplings 


K = iJt, Ki = iJj_t, (6) 


which in the present nonequilibrium context, however, 
are now complex. Notice that the above identification 
of Q(t) with a partition function is independent of di¬ 
mension. In the following, the ID and 2D cases will be 
considered because they allow for exact solutions. 


One dimension:- The ID Ising model can be solved on 
the basis of the transfer matrix T : (| 

6«) = «rT-, r=‘( e i ';**), (7) 

with tr denoting the trace over a basis of the 2x2 matrix 
problem. The matrix T has two eigenvalues v c = ch(A) 
and v s = sli(A') with sh(A') and ch(A) the hyperbolic 
sine and cosine, respectively. For L —> oo, Q(t) is 
dominated by the eigenvalue of largest magnitude, i.e., 
Q(t) = v L with v = v c if \v c \ > \v s \ and v = v s oth¬ 
erwise. Notice, that v can switch between v s and v c 
yielding a nonanalytic structure in A (t) = — 2Re [log(^)]. 
This switching of the dominant eigenvalue is the under¬ 
lying origin of the DQPTs in ID Ising models, as has 
also been seen in XXZ chains [71. The critical times t n 
of the DQPTs are given by the condition \v c \ = \u a \, i.e., 
t n = 7t(2 n + 1)/(4J) with n € Z, compare also Ref. [2]. 
In equilibrium, \v c \ = \v s \ can only be satisfied in the 
limit of zero temperature T = 0. When discussing the 
anticipated RG procedure below it will be shown that 
this correspondence is not accidental but rather has a 
very profound origin. 

Having established the presence of DQPTs it is now 
the aim to address the question of scaling and universal¬ 
ity. For that purpose an RG scheme in complex param¬ 
eter space will now be introduced. RG transformations 
in complex parameter space have been previously stud¬ 
ied in the context of the standard model era as well 
as for equilibrium partition functions in complex param¬ 
eter spaces m 122j . Eliminating every second spin via 
decimation [271(23] an exact RG transformation can be 
formulated yielding the following recursion relation [24] 

th(A') = th 2 (A). (8) 

As a result, it is found that the RG has two fixed points 
K* = 0, oo corresponding to the equilibrium ones at 
infinite and zero temperature even when K is complex 
initially. For K with A' <C 1 this gives K' = K 2 
implying that the fixed point A* = 0 is stable. For 
A = K* + SK in the vicinity of K* = oo one obtains 
that SK' = 2 SK = b x SK with 6 = 2 the change in 
length scale due to the decimation and A = 1 the asso¬ 
ciated anomalous dimension. Thus, it is found here that 
the fixed point K* = oo is unstable as in the equilib¬ 
rium case. But remarkably, this is not necessarily true 
for initial couplings beyond the linear regime. In partic¬ 
ular, the DQPTs at times t n map onto the K* = oo fixed 
point after precisely two RG steps. Times t with weak 
deviation r = (t — t c )/t c from a DQPT at t c map after 
two RG steps onto the linear regime of the unstable fixed 
point. Using that A = 1, one can then directly deduce 
the scaling form of g(t): 

d( T ) ~ |r| d/A 4>± = |t|$±, T= t —^. (9) 
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with dimension d = 1 and a constant that may differ 
for r ^ 0. Remarkably, this is indeed the scaling behavior 
that is found from the exact solution of g(t) in the vicin¬ 
ity of the DQPTs [2], compare also Fig. [lj Therefore, the 
nonanalytic behavior of g{t) can now be attributed to an 
unstable fixed point of an RG allowing to extend funda¬ 
mental concepts such as robustness, scaling, and univer¬ 
sality to the nonequilibrium regime on general grounds. 
Notice that robustness has been established for particular 
cases recently mm- In the present work a specific initial 
state and final Hamiltonian has been considered so far. 
However, the identification of DQPTs with unstable fixed 
points allows to conclude that weak symmetry-preserving 
perturbations do not change the universal properties. In 
the following this will now be demonstrated by incorpo¬ 
rating a weak transverse field in the final Hamiltonian. 

Transverse Fields. - For h/J <C 1 the field part V = 
~h Zi°f of Eq. § can be eliminated using standard 
time-dependent perturbation theory [24]. To first or¬ 
der in h/J one obtains that Q(t) = 2~ L Tre H can be 
again represented in terms of an effective classical Ising 
model H , but now including also next-to-nearest neigh¬ 
bor (NNN) interactions [24] : 

H = KY J «+i + G Y, * 1 °!+ 2 ( 10 ) 

i i 


with G = —iht/2 + ihsm(4Jt)/(8J) and a modified cou¬ 
pling K = iJt + h[ 1 — cos(4Jt)]/(4J). Within the same 
decimation RG as used for the NN case let us eliminate 
every second lattice site. Based on a cumulant expan¬ 
sion [25] for the perturbative NNN couplings one obtains 
to first order in G: [24] 


K' = P + G 


i —4 pi 

1 — e 

+ 


G' = G 



( 11 ) 



Figure 2: (color online) Dynamics of A (t) in the 2D Ising 
model, (a) For strong anisotropies j± = J±/J -C 1 the weak 
coupling J_l represents an irrelevant perturbation and the crit¬ 
ical properties, i.e., the kinks, are identical to the ID case, 
(b) For the isotropic 2D Ising model the nonanalytic struc¬ 
ture changes to a logarithmic singularity which is illustrated 
in the second derivative shown in (c). 


with e q the solution of the equation: 

ch (e q ) = ch(2A'_L)ch(2/i) — sh(2A'j_)sh(2AT) cos(q). (13) 

and s(x) = sgn[7£(x)] returns the sign of its argument’s 
real part. Here, K is given by the condition th(7F) = 
exp(—2 K) with th(AT) = sh(A')/ch(AT). 

In Fig. [2j the dynamics of A (t) is shown for different 
anisotropies j± = J±/J. For jj_ <C 1, DQPTs are found 
at times t n = (2 n + l)7r/(4J) which are solely controlled 
by the coupling J. Indeed, it will be shown below using 
a perturbative RG that a weak coupling j± <C 1 repre¬ 
sents an irrelevant perturbation. If, however, j± = 1 a 
drastic change in the nature of the DQPT occurs with a 
logarithmic nonanalyticity: 

g(r) ~ t 2 log(|r|), (14) 


with th(P) = th 2 (A') the solution at h = 0, see 
Eq. ([8|. This set of RG equations exhibits two fixed 
points ( K*,G*) = (0,0), ( 00 , 0). In the vicinity of the 
unstable fixed point K* = 00 we get SK' = 2SK + 3G/2 
and G' = G/4 and therefore weak fields ft/J <C 1 consti¬ 
tute an irrelevant perturbation. This is in perfect agree¬ 
ment with the exact solution where it can be seen that 
h/J > 1 is necessary to destroy the DQPTs [2j, see also 
Fig- 0 Moreover, the scaling properties of the DQPTs 
are invariant under a slight modification of the initial 
state by taking the ground state for an initial transverse 
field 1 < h/J < 00 . Therefore, it is expected that the 
main results are also valid beyond the case of a fully po¬ 
larized state studied here. 

Two dimensions:- The partition function of the 2D 
Ising model can be solved exactly [2141251 . For the case 
of complex AT this is still possible which yields: 

g(t) = - * log[2sinh(/0] - J ^s(GjH (12) 


which is illustrated in Fig. [2j Notice the remarkable simi¬ 
larity of the scaling behavior in Eq. (14) with the equilib¬ 


rium free energy at the equilibrium critical point of the 
2D Ising model when r denotes the relative temperature 
distance to the critical point [2T|. 

As opposed to ID, it is not possible to derive a closed 
set of exact RG recursion relations for the 2D case. In the 
limit of strong anisotropy j± <C 1, however, an approxi¬ 
mate RG can be constructed. For that purpose let us de¬ 
compose the square lattice into even and odd rows. The 
odd rows can be eliminated perturbatively using a cumu¬ 
lant expansion 125] which in the present case is controlled 
via \K±\ -C 1. To second order in K± one obtains: [24] 


K' = K + 2QK] , K\ = K\ 


(15) 


Here, Q = th(It) if \v c \ > \v s \ and Q = l/th(A') other¬ 
wise with v c and v s the eigenvalues of the ID T-matrix, 
see Eq. 0- For A'_l < 1 initially, K* ± = 0 is always ap¬ 
proached implying that K± is an irrelevant perturbation. 
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Figure 3: (color online) (a) Dynamics of C z (t) for h/J = 0.1 
on a N = 5 x 5 square lattice for different anisotropies j±. 
(b) Power-law scaling of x(Y) in the vicinity of the DQPT. 
Brown curves are on the NN square lattice for j± = 1/8 
with d* = 1 (upper curve) and the isotropic limit j± = 1 
with d* = 2 (lower curve). The dashed blue curves for the 
long-range Ising models relevant for trapped ions in ID (up¬ 
per curve) and in 2D (lower curve) are also included. As a 
reference, quadratic and quartic power-laws are also shown 
demonstrating the dynamical scaling of Eq. (171. 


As a consequence, the fixed point describes a set of un¬ 
coupled ID chains. Indeed, \(t) displays kinks, confirm 
Fig- [2} In this context let us therefore introduce an ef¬ 
fective dimension d* which takes a value d* = 1 for the 
ID system and also for the strongly anisotropic 2D one. 

Decreasing the anisotropy makes the RG transforma¬ 


tion in Eq. (151 less and less controlled. In particular, the 


isotropic point J± = J and its associated logarithmic sin¬ 


gularity of Eq. (14) is not accessible in this way. Thus, 
within the current methodology, a rigorous identification 
of this DQPT with a fixed point is not possible. The par¬ 


ticular scaling form of g(t) in Eq. (14), however, suggests 
that the DQPT in the isotropic limit is controlled by the 
unstable fixed point of the 2D Ising model. A further ar¬ 
gument supporting this hypothesis is given below when 
discussing the power-law scaling of the spin correlations. 

Spin correlations:- Having unraveled scaling and uni¬ 
versality of the Loschmidt amplitude, the major chal¬ 
lenge now is to connect to local observables. Although for 
particular quenches out of symmetry-broken phases such 
a connection has been established a general under¬ 
standing, however, is still lacking. As the present DQPTs 
are associated with the equilibrium critical points of Ising 
models, it is a question of fundamental importance of 
how the divergent equilibrium correlation length becomes 
manifest in the nonequilibrium dynamics. Due to Lieb- 
Robinson bounds it is, of course, not possible to build 
up diverging correlations within a finite time interval. 
Nevertheless, it will now be demonstrated that the un¬ 
derlying DQPTs are responsible for the buildup of NN 
spin correlations within the rows: 


Im 


(16) 


In Fig. [3] the dynamics of C z (t) is shown. Notice that 
the case h = 0 constitutes a singular limit because there 
C z (t) is a constant of motion such that C z (t ) oc h/J 
for h/J -C 1 and C z (t)J/h becomes a universal function 
independent of h in the limit h —> 0. As one can see, 
the spin correlations develop a maximum in the vicinity 
of the DQPT at t c = 7t/(4J) with a slope, however, that 
differs for the isotropic 2D case compared to those with 
effective dimension d* = 1. This is associated with a 
different power-law scaling that can be quantified via: 


x(t) = ^ [C z (t c ) - C z {t c + r )] h cx° r 2d \ (17) 


see Fig. [3] This scaling only depends on the effective 
dimension d* of the DQPT and therefore only on the 
DPQT’s universality class if we assign d* = 2 for the 
DQPT satisfying the scaling of Eq. (14) for g{t) equiv¬ 
alent to the critical point of the 2D Ising model. The 
numerical data in Fig. [3] has been obtained from exact 
diagonalization (ED) using a Lanczos algorithm with full 
reorthogonalization EDI. 


Trapped ions:- This dynamical scaling can be observed 
experimentally in systems of trapped ions within the cur¬ 
rent technology as will be outlined in the following. Fully 
polarized states as required in Eq. Q can be initial¬ 
ized with a high fidelity m ■ Coherent time evolution 
of transverse-field Ising Hamiltonians has been demon¬ 
strated both for ID [5TH55] and 2D|M]. The Ising cou¬ 
plings, however, are not of NN type but rather long- 
ranged, Jim = J/\n — r m \ a with 0 < a < 3 [35: and 
r; the location of the ion in real space. In Fig. [3j nu¬ 
merical ED data of y(r) for the long-range potentials is 
included for the dipolar case a = 3 with lattice spacing 
a = 1. Here, open boundary conditions have been used 
and x(r) includes only those NN correlations which do 
not contain spins at the boundary to minimize boundary 
effects. As the simulations indicate, the spin correlations 
also obey the dynamical scaling of Eq. making it 
accessible within current trapped ion technology. 


Conclusions:- It has been shown that DQPTs in ID 
and 2D Ising models are controlled by unstable fixed 
points of complex RG transformations opening the pos¬ 
sibility to apply the concepts of scaling and universal¬ 
ity to the out-of-equilibrium regime. Importantly, this 
leads to a dynamical scaling of the spin correlations which 
only depends on the universality class of the underlying 
DQPT and which is accessible experimentally in systems 
of trapped ions. 
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Supplemental Material to 

Dynamical quantum phase transitions as critical points of complex 
renormalization group transformations 


In this supplemental material, we provide methodological details for the derivation of the 
real-space RG equations. 


A. REAL-SPACE RG IN ID: DECIMATION 

In the following, the well-known real-space decimation scheme for the one-dimensional Ising model is recapitulated, 
see Refs. flTWIli . This will be the basis for the analysis below. Consider the (complex) partition function Z(K,N) of 
the one-dimensional Ising model: 


Z(K, N ) = Tr 


e H(K,N) 


N 

H(K,N)=kJ2<7?<7i+i 

1=1 


(1) 


of length N with K a complex coupling in general, see main text. For simplicity, periodic boundary conditions, 
i.e., (J^r +1 = erf, are used. Within decimation, it is the aim to eliminate every second spin yielding Z(K,N) = 
Z(K ', N/ 2) = Trexp[H(K', N/ 2)] with a new Hamiltonian H(K ', N/ 2) for the remaining N/2 spins that now interact 
via the renormalized coupling K'. As the Pauli matrices erf mutually commute, the partition function can be factorized 
such that integrating out all the odd spins gives: 

N N/2 _ 

Z(K,N) = Tr e Tr 0 JJ e K<7 ‘ v ‘+ 1 = Tr e ^ e^> 2i+1 (2) 

1=1 1=1 CT~ l + 1 =±l 


with Tr e / 0 denoting the trace over the even (e) or odd (o) sites, respectively. Using exp [Kaa 1 ] = cosh(A') + sinh (K)cra' 
for cr and er' Pauli matrices, one obtains: 


y e Kcr 2i^2i+i e K^ l+1 cr^ +2 _ cos ]^^x) [l + tanh 2 (/i')crf;crf; +2 ] , 

a 2l + l = ±l 

This gives the desired RG equation for the couplings: 

tanh(A'') = tanh 2 (AT), e 2K = cosh(2 K). 


(3) 

(4) 


Having recapitulated the decimation RG scheme for the classical ID Ising chain, it is the aim of the following to study 
the influence of weak transverse fields. 


B. WEAK TRANSVERSE FIELDS 

In this section the perturbative calculations for the Ising chain in presence of transverse fields are shown. First, 
the mapping onto the classical Ising with next-to-nearest neighbor interactions is presented using time-dependent 
perturbation theory. Afterwards, the RG equations of the main text are derived on the basis of standard cumulant 
expansion. 


Mapping onto classical Ising chain 

The Ising chain with a transverse field h reads: 

H = H 0 + V, H 0 = ~jJ2 V=-hY,°i- 

l l 


( 5 ) 
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In order to determine the Loschmidt amplitude Q(t) = (^o| exp[—see Eq. (1) of the main text, for a weak 
transverse field, let us turn to an interaction picture with respect to Hq : 

e~ iHt = e~ iHot W(t), W(t) =Te~ i ti dt ' v{t '\ V{t) = e iHot Ve~ iHot , (6) 

with T the usual time-ordering prescription. For the Loschmdit amplitude this implies: 

G{t) = {-i/j 0 \e~ lHoi W(t)\-i/) 0 ) =^2e~ lEst (i/j 0 \s)(s\W{t)\ip 0 ), (s\W(t)\ip 0 ) = {s\Te~^o |^ 0 ). (7) 


by inserting an identity via 1 = |s)(s| with |s) = |si,... ,sl) denoting the complete basis of spin configurations 

and si = dzl the orientation of the local spin at site l. Using standard cumulant expansion methods for time-ordered 
exponentials m one obtains to lowest order in the transverse-field strength h: 


(s |Te l ti dt \ip 0 ) « (s|V’o)exp 
For the choice |^ 0 ) as the ground state of V of the main text this gives 
(s\Te~ l f° dt \ip 0 ) = 2 -i/2 exp[/i(l - cos(4 Jt ))/(4 J) 7 sis i+1 

i 



(s\v(mo y 

(s|^o) . 


- - sin(4Jt)/(4J)) 7 s z s ;+2 ] 

l 


( 8 ) 


(9) 


This result can be rewritten in terms of a classical Ising model but now including also next-to-nearest neighbor (NNN) 
interactions which in the end gives for the Loschmidt amplitude the desired result: 


g(i) = Tre ff , H = KJ2 ^X+i + G 7 of +2) 

l l 


with 


K = iJt + 77 [1 — cos(4 Jt )], G = — i^-t + sin(4Jt). 

4 J 2 8 J 


( 10 ) 


( 11 ) 


RG equations 

It is the aim of the following to apply the real-space decimation approach to the Ising chain with NNN interactions. 
Again, we will eliminate every second spin, which we take to be odd lattice sites without loss of generality: 


.H 


= Tr 0 e ff , H = H 0 + V, H 0 = af<rf+i* V = G 7 < 


l a l+ 2- 


( 12 ) 


Here, H denotes the renormalized Hamiltonian after this RG step for the remaining even lattice sites and Tr 0 is the 
trace over all the odd lattice sites that are supposed to be integrated out. As [Hq, V] = 0 it holds that 


e H = Tiv, 


.H 0 


Tr„ 


e Ho e V 


Tr 0 e H ° 


= Ttv, 


.Ho 




(13) 


with ((...)) denoting a generalized average. As V is proportional to G, which is assumed to be small, one can perform 
a cumulant expansion |25] which to lowest order is 


((e v )) « exp [((V))] 


( 14 ) 


Evaluating (( V )) yields the desired RG equations presented in the main text. 














C. ANISOTROPIC 2D ISING MODEL 


The anisotropic 2D Ising model considered in the main text reads: 

H = a l,m a l,m +1 — J-L a l,m a l+l,m 

lm lm 


(15) 


with l ( m ) denoting the row (column) on the 2D square lattice. The coupling within the rows is of strength J whereas 
along the columns it is J±. In the following, a perturbative real-space RG will be derived taking j± = J±/J as a 
small parameter. Each isolated row is an exactly solvable ID Ising model. It is the aim the eliminate every second 
(odd) row. For the Loschmidt amplitude Q(t) one can then introduce the decomposisiton 

6(t) = tv [e H ], n = n e + n 0 + 'H e0 (16) 

with 


H e = Kj2 l 


H-o — K C^Z+l.mO^+l.TO+l' H eo — K <J l,m a l+l.m- 


(17) 


lm lm Ir, 

The elimination of the odd rows can be formally implemented in the following way: 


, Tr \e H o e H co ] 

e n = e He Tr 0 [e w »] ^ ^ J = e«=Tr 0 [e«»] «e w “)). 

Because T-Leo is proportional to the weak coupling K± -C 1 it is again possible to perform a cumulant expansion: 


(18) 


«e*“»«exp 


2 {{HeoHeo}) 


(19) 


which now has to be determined up to second order because the first order vanishes. Neglecting the perturbatively 
generated longer-ranged interactions one obtains the desired RG equations 


K' = I< + 2 QKl , K ' L = K \, 


( 20 ) 


with Q the nearest-neighbor correlation function of the ID Ising model which is Q = tanh(A') if \i/ c \ > \v a \ and 
Q = l/tanh(A') otherwise [TJ. 





